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We propose a technique to realize a tachyonic band structure in a coherent network, such as an 
array of coupled ring resonators. This is achieved by adding “PT symmetric” spatially-balanced gain 
and loss to each node of the network. In a square-lattice network, the quasi-energy bandstructure 
exhibits a tachyonic dispersion relation, centered at either the center or corner of the Brillouin zone. 

There is one tachyonic hyperboloid in each gap, unlike in PT-symmetric tight-binding honeycomb 
lattices where the hyperboloids occur in pairs. The dispersion relation can be probed by measuring 
the peaks in transmission across a finite network as the gain/loss parameter is varied. 


Most wave theories, including but not limited to 
quantum mechanics and classical electromagnetism, are 
formulated using equations of motion with Hermitian 
Hamiltonians. In quantum mechanics, Hermiticity en¬ 
sures the general conservation of total probability under 
time evolution; conversely, its violation describes am¬ 
plification (gain) and/or loss. Thus, for instance, non- 
Hermitian Hamiltonians are used in effective theories of 
decaying quantum systems, in which the wavefunction 
can leak away into unmonitored degrees of freedom. In 
optical physics, gain and loss processes are even more 
ubiquitous, in the context of the emission and absorption 
of light, and the prescriptions for dealing with these pro¬ 
cesses (e.g., introducing complex frequency-domain di- 
electic permittivities) are similarly well known. 

Several years ago. Bender and co-workers made 
the striking observation that in systems possessing 
parity-time {PT) symmetry, corresponding to spatially- 
balanced gain and loss, the Hamiltonian can have purely 
real eigenvalues (i.e., probability-conserving eigenstates) 
despite being non-Hermitian mm- Subsequently, a se¬ 
ries of works showed both theoretically and experimen¬ 
tally that this effect could be demonstrated in optical 
structures, using optical gain and loss [3HS1 • In PT sym¬ 
metric optical lattices lilSl, the photonic band structure 
has quite unusual features: the band energies can be real 
in one region of the Brillouin zone, where the Bloch eigen¬ 
states are PT symmetric, and complex in another region 
where the PT symmetry is spontaneously broken. 

For 2D lattices, Szameit et al. showed that the PT 
symmetry-breaking phenomenon has a startling interpre¬ 
tation in terms of emergent “tachyons”: hypothetical su¬ 
perluminal particles which are not known to exist in na¬ 
ture [30]. A two-dimensional honeycomb lattice can be 
realized using an array of coupled optical waveguides. In 
the Hermitian case, the bandstructure is graphene-like, 
featuring a pair of linear band crossing points (“Dirac 
points”) with band velocity vd [inilll]- When gain and 
loss are added to alternating sites of the honeycomb lat- 
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tice, the Bloch Hamiltonian becomes non-Hermitian, and 
in the vicinity of each Dirac point it takes the form of a 
Dirac Hamiltonian with imaginary mass. The eigenstates 
are tachyons whose group velocities are larger than vd- 
In fact, the group velocities become infinite along a “crit¬ 
ical” ring in fc-space surrounding each Dirac point, corre¬ 
sponding to the PT symmetry breaking transition points 
of the Bloch Hamiltonian. A similar dispersion relation 
has also recently been realized in Ref. m- However, 
it should be noted that the notion of group velocity as 
the slope of the dispersion relation must be reevaluated 
in non-Hermitian systems. Using the Hellman-Feynman 
theorem, it has been shown that as the critical ring is 
reached in k-space, significant corrections to this defini¬ 
tion of the group velocity arise [32H34] . 

This paper describes an alternative way to realize a 
tachyon-like bandstructure, using a 2D network of coher¬ 
ent waves [HIS] with non-unitary evolution. A “net¬ 
work model”, unlike the tight-binding models commonly 
used in condensed-matter and optical physics, does not 
describe a lattice in terms of a Hamiltonian. Instead, 
it uses an evolution matrix to describe the propagation 
of waves through a network of directed links and nodes. 
As discussed below, such networks can be realized in a 
variety of ways, such as coupled optical resonator lat¬ 
tices [IM[, microwave networks [H [21] , and RF cir¬ 
cuits [2S]. Network models can produce bandstructures 
with various unusual features that are not found in static 
Hamiltonian models |15j ; in fact they can be mapped to 
the class of “Floquet” systems, described by Hamiltoni¬ 
ans that vary periodically in time [TTII^HUSH] . As we shall 
see, introducing PT symmetric gain and loss to a square- 
lattice (not honeycomb) network yields a bandstructure 
with tachyonic Dirac dispersion relations. But unlike 
the previously-studied tight-binding honeycomb lattice, 
where the tachyonic Dirac hyperboloids occur in pairs, 
this network bandstructure contains a single hyperboloid 
in each gap. Finally, we will show how the tachyonic 
dispersion relation’s critical /c-vector can be measured 
through transmission experiments across finite networks. 

Consider the network model shown schematically in 
Fig.g It consists of links and nodes, where each link car- 
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FIG. 1. (a) Schematic of a chiral square-lattice network. 

Waves propagate directionally in the links (arrows), form¬ 
ing clockwise loops, and couple at the nodes where adjacent 
loops approach each other (shaded ovals). Two of the loops 
are highlighted in red and blue for ease of comparison with 
(b) and (c). (b) Detail of a unit cell of the network, showing 
the wave amplitudes entering and leaving each node. These 
amplitudes enter into Eq. 0- 


ries a one-directional wave described by a complex scalar 
amplitude; these are arranged in a 2D square lattice, with 
each cell containing four links arranged in a chiral loop 
[T^ . Adjacent loops are coupled at the nodes, which are 
described by 2 x 2 scattering matrices. This model was 
first introduced for studying the transport properties of 
disordered quantum Hall systems m-, it captures the es¬ 
sential features of a disordered 2D electron gas in strong 
magnetic fields, where the electronic orbits follow chiral 
“race-tracks” along the equipotentials of a disordered po¬ 
tential landscape, and can tunnel to adjacent race-tracks 
at potential saddle-points [12]. There is now an exten¬ 
sive literature on the use of network models for studying 
electronic transport; see Ref. [13] for a survey. 

Recently, researchers have implemented chiral net¬ 
works in classical electromagnetic settings. One type of 
realization is an on-chip coupled resonator lattice [TOlllSj . 
of the sort proposed and experimentally studied by Hafezi 
et al. Optical ring resonators are arranged in a lattice, 
playing the role of the network’s chiral loops. Each pair 
of adjacent resonators is coupled by an auxiliary ring 
waveguide, which acts as a node. Due to local momentum 
conservation at the inter-waveguide interfaces, the opti¬ 
cal modes of the lattice decouple into one set of modes 
where light propagates clockwise in the main rings, and 
another that is counter-clockwise; each set maps onto a 
network model. Such resonator lattices can exhibit topo¬ 
logical edge states and fractal Hoftadter spectra [Ta[iH], 
as well as topological transitions and anomalous topo¬ 
logical phases [TllEIlEaEl]. A chiral network can also 
be realized using a microwave circuit |23j . The nodes of 
the network are implemented using directional couplers; 
auxiliary rings are not necessary, since the microwave 
components need not be strictly planar. The chirality of 
the network can be enforced using microwave isolators. 

Regardless of the network model’s underlying imple¬ 
mentation, its properties can be described theoretically 
in terms of evolution matrices. And for a disorder-free, 


spatially infinite network, the evolution matrix descrip¬ 
tion gives rise to a bandstructure [H dg. We briefly 
review the procedure. Fig.j^b) shows a schematic of one 
cell of the network, which contains two couplers connect¬ 
ing adjacent loops along the x and y directions. These 
couplers can be described by scattering matrices Sx and 
Sy. We denote the four input wave amplitudes into these 
couplers by { 6 i,..., 64 }, and the wave amplitudes on the 
other side of those links by {oi,..., 04 }. We assume each 
link has equal phase delay </>, so that bn = For the 

moment, we ignore gain and loss, so that Sx and Sy are 
unitary and (/) is real. Using Bloch’s theorem, we can 
relate the input and output amplitudes as follows: 


Sx 

'hi 

.^2. 

= 

036-*'=-' 

046®*’“ 


Sy 

^3 

= 

0.2 

ky) 


where k = [kx, ky] is the Bloch wavevector, with the lat¬ 
tice spacings normalized to 1. These can be combined 
into a single 4x4 eigenvalue equation: 
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For each k, U{k) is unitary, and its eigenvalues are 
exp(—z^), where the “quasi-energies” (j) are the discrete 
values of the link delay for which modes can propagate 
with the given k. Because the bandstructure is defined by 
an evolution operator, U{k), rather than a Hamiltonian, 
it falls into the same class as the “Floquet” bandstruc- 
tures describing periodically-driven lattices mmm- 
Note that is an angle variable, unlike the energy occur¬ 
ring in a conventional Hamiltonian eigenproblem. 

The bandstructure depends on the choice of the 2x2 
unitary matrices Sx and Sy, each parameterized by four 
Euler angles. It turns out that the bandstructure topol¬ 
ogy is determined only by one pair of Euler angles, de¬ 
noted by Ox^y in Ref. m, which parameterize the cou¬ 
pling strengths between adjacent loops. The gaps in the 
quasi-energy bandstructure close when 9x + 0y = 7!'/2 
[2T] . For simplicity, we fix the other Euler angles so that 




cos Ofi i sin 0^ 
i sin 9y cos 0^ ’ 
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This describes a coupler which is symmetric under 180° 
rotations [see Fig. [^b)], and behaves the same when the 
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FIG. 2. Bandstructure of the network model for 9^ = By = 
7r/4. There are Dirac points at 0 = {—7r/4, 37r/4} and = 
ky = 0, and at 0 = {—37r/4,7r/4} and k^ = ky = ±7r. 


order of the two inputs, and the two outputs, are simul¬ 
taneously swapped. 

Fig. shows the spectrum for dj; = 9y = njA. There 
are four bands, joined by Dirac points at 

( ka: = ky =0, 

\ G {-7r/4,37r/4} 

f kx — ky = TT, 

[ (/) G {—37r/4, —7r/4}. 


These Dirac points can be conveniently derived by taking 
the squared matrix mi: 


U\k) = 


'Ui{k) 0 
0 U2{k) 

Ur = WySyWxSy,, 


( 8 ) 


U2 = Wy,Sy,WySy. 


C/^ consists of two blocks with identical spectra, and its 
eigenvalues are exp(—210 -|- 2rmr) for m G Z. Focusing 
on the first block, we let = tt/ 4 -b A9, and expand to 
first order in ky, and A9. The result is 


« exp i Hoik) 


TT 

2 J 




}, (9) 


where 


HD{k) = '^^^^kx,+ '^^L^ky-2A0ax. (10) 


This is a Dirac Hamiltonian with mass 2A0 and band 
velocity vd = 2“^/^. Another pair of Dirac points is 
obtained by expanding ky = tt + Ky. The result is 


Ui 


exp < —i 


HD{n) + - 



( 11 ) 


which gives the second set of Dirac points in Eq. Q. 

In the critical bandstructure {Ad = 0), there is a sin¬ 
gle Dirac cone in each gap. This contrasts with the 
more familiar case of time-reversal symmetric honeycomb 


lattices (e.g. graphene), where the Dirac points occur 
in pairs, which is a manifestation of the “fermion dou¬ 
bling” principle |29j . The unpaired Dirac point in the 
network model’s bandstructure is reminiscent of the un¬ 
paired massless chiral relativistic fermions which occur 
at the Haldane model’s critical points under broken in¬ 
version and time-reversal symmetry |29j . In the network 
model, those symmetries are likewise explicitly broken. 

Interestingly, this critical bandstructure cannot be de¬ 
scribed adequately by an effective Hamiltonian Heg{k) = 
jlog[17(A;)]. This is because the quasi-energy bandstruc¬ 
ture is completely ungapped; the eigenvalues of C/(fc), 
over all k, cover the unit circle. Thus, there is no way to 
assign the logarithm’s branch cut, without there being a 
locus of ^-points where some eigenvalues of U{k) cross 
the cut. There is no Hesik) which can be “smoothly” 
defined over all k in the Brillouin zone. 

If we tune the coupling strengths away from the critical 
point, the quasi-energy bandstructure becomes gapped. 
For A9 < 0, it is topologically trivial, and for A0 > 0 it 
is topologically non-trivial [HI [21]. In both cases, how¬ 
ever, it can be shown that every band has zero Chern 
number m- This happens because, in the critical band- 
structure, each band had a Dirac point above and below, 
a situation that is possible because 0 is an angular vari¬ 
able and hence not bounded above or below. Hence, for 
A9 > 0 the network is in an “anomalous Floquet insu¬ 
lator” phase, exhibiting topological edge states despite 
all bands having zero Chern numbers. Similar anoma¬ 
lous phases are also known to occur in periodically-driven 
Floquet topological insulators [HHISj. 

We are now ready to consider a network contain¬ 
ing gain and loss. In the honeycomb lattice, Szameit 
et al. have previously shown that adding PT symmet¬ 
ric gain and loss to the alternate sublattices distorts 
the bandstructure’s Dirac cones into hyperboloids, corre¬ 
sponding to two species of emergent tachyons [30|. This 
is caused by the Bloch states near each Dirac point un¬ 
dergoing spontaneous PT symmetry breaking. In the PT 
symmetric region, the bands are real and have group ve¬ 
locity exceeding the Dirac velocity vjj. The group veloc¬ 
ity approaches infinity at the waists of the hyperboloids, 
which are the PT symmetry breaking points of the Bloch 
Hamiltonian. 

In the network model, tachyonic behavior can arise by 
setting A9 to be imaginary. This can potentially be real¬ 
ized in the optical resonator domain by using ‘auxiliary 
rings’ that lie in between the principal rings, and are opti¬ 
cally pumped and thus have gain (unpumped rings would 
naturally exhibit loss). A candidate platform would be 
that used in Ref. m- In the context of microwave net¬ 
works, auxiliary directional couplers could be used in a 
similar way, in combination with amplifiers. According 
to Eq. (101, this gives the effective Dirac Hamiltonian 


an imaginary mass. For A9 = ij, the coupling matrices 
become 
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FIG. 3. Section of the network bandstructure for 0^ = 6y = 
(0.25 + 0 . 02 i) 7 r. Each of the Dirac points becomes a hyper¬ 
boloid. We have zoomed in on the hyperboloid at ~ — 7 r/ 4 , 
^ ^ 0 . 


This yields the bandstructure shown in Fig. Each 
Dirac cone becomes a hyperboloid, corresponding to a 
tachyonic dispersion relation. Since there was originally 
only one Dirac cone per gap, the hyperboloids are un¬ 
paired, unlike in the PT symmetric honeycomb lattice 
m- The band quasi-energies are all real, except for the 
regions of k inside the waists of the hyperboloids. Using 
Eq. (10), we find the critical wavenumbers 


kcij) = 2^/^7 -b 0(7^). 


(13) 


The coupling matrix of Eq. ( |l2| ) has the same 180° 
rotational symmetry as the previously-discussed unitary 
coupling matrix of Eq. ([^ . However, for 7 0 it is man¬ 

ifestly non-unitary. This may be seen from the eigenval¬ 
ues a± = a ± ia*, whose magnitudes are |cr±| = e^'^. 
The corresponding eigenvectors are [1;±1]; one of these 
eigenvectors is amplified, and the other is damped by an 
equal and opposite amount. This is very similar to the 
behavior of scattering matrices derived from the wave 
equation in PT symmetric media |351 - I37] . Furthermore, 
the coupling matrix can be decomposed as 


S^.=So 


( 7 + 0 

0 (T_ 


So, where Sq = 


1 1 

\/2 

V2 


■ (14) 


Thus, such a coupler can be implemented by passing the 
inputs through a unitary 50:50 coupler described by 5o, 
applying balanced gain and loss to the results, and then 
re-mixing through a second Sq coupler. 





—TT —7r/2 0 n/2 tt —7t/2 0 7r/2 

Re(fc3,) Im(fcj,) 


FIG. 4. (a) Schematic of a semi-infinite network “strip” of 

width N in the y direction, (b)-(c) Band diagram of real 
(p versus complex kj,, with periodic boundary conditions ap¬ 
plied at the upper and lower edges of the strip. Here, we 
take A'’ = 10 and coupling parameter 9 = 7 r /4 -|- O.li. At 
4> = — 7 r/ 4 , the red arrow indicates a pair of tachyon modes, 
which have Re[fca;] 7 ^ 0 (left plot) and Im)^^;] = 0 (right plot). 
At the same quasi-energy, there is a set of evanescent modes 
with Re[fca;] = 0 and Im[fe 2 :] yf 0. (d)-(e) Band diagram with 
Dirichlet boundary conditions at the strip edges. The Dirich- 
let boundary conditions spoil the PT symmetry, causing the 
tachyonic states to have small but non-zero Im[A:i,]. 


How might the tachyonic bandstructure be experimen¬ 
tally verified, whether in the context of optical ring res¬ 
onators or microwave networks? One possibility is to 
construct a wavepacket and show that its group ve¬ 
locity can exceed the effective Dirac velocity, as dis¬ 
cussed in Ref. [SD]. For the network model, this ap¬ 
proach could work if (j) is proportional to frequency, and 
the coupling parameters are approximately frequency- 
independent |22j . However, the natural quantities to 
study in a network model are the steady-state reflection 
and transmission for fixed tp. Here, we present an al¬ 
ternative experimental approach for probing the tachy¬ 
onic bandstructure, based on measuring the transmission 
across a set of finite networks. 

Let us first examine the bandstructure of a “strip” of 
network, shown schematically in Fig.|^a). The strip ex¬ 
tends infinitely in the x direction, and has a width of 
N cells in the y direction. There are two useful choices 
of boundary conditions that we can impose on the edges 
of the strip. Firstly, we can impose periodic edges by 
making row -|- 1 equivalent to row 1 (i.e., rolling the 
strip into the surface of cylinder). Secondly, we can im¬ 
pose “Dirichlet” edges by terminating the network at the 
edges of the strip, setting = 04 and ■ (We 

could also introduce phase factors into these edge rela¬ 
tions; but that generates additional non-topological edge 
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states, which we are not interested in here.) 

Fig. ID shows the band diagram of real-()) versus 
complex-fca;, for the semi-infinite strip. This band dia¬ 
gram is calculated from the eigenvalues of the transfer 
matrix across one x period of the strip m- The rea¬ 
son for plotting real (j) versus complex kx, rather than 
complex (j) versus real kx, is that we will be interested 
in the propagation of modes at a fixed real quasi-energy 
(p, chosen to correspond to one of the Dirac points. We 
first focus on Fig. |^b)-(c), which shows the case of pe¬ 
riodic edges. For each real p, all the modes are either 
purely propagating (real kx), or purely evanescent (imag¬ 
inary kx)- The evanescent modes are completely non¬ 
propagating (Re[fca;] = 0), and are similar to the evanes¬ 
cent modes which occur within the band gaps of ordinary 
Hermitian systems. As for the propagating modes, there 
are specific branches of these modes which have tachyonic 
dispersion relations, and are highlighted in red and gold 
in the figure. As indicated in Fig.j^c), these modes prop¬ 
agate with no amplification nor dissipation (Im[fca;] = 0). 

Fig. gd) -(e) shows the band diagram for Dirichlet 
edges. In this case, the modes are no longer purely prop¬ 
agating or purely evanescent, but have complex kx- This 
happens because the edge conditions spoil the PT sym¬ 
metry of the network. Nonetheless, the projected band 
diagram remains qualitatively similar to Fig. ib)-(c). 
In particular, there are tachyon modes which are weakly 
damped (small |Im(fca;)|) compared to the other modes. 

We can now formulate a transmission-based experi¬ 
mental signature for the tachyon modes. Consider a strip 
of length M (in the x direction). Along column I, we in¬ 
ject equal wave amplitudes a" = into each of the 

rightward links [see the schematic in Fig. |^a)]; then, M 
columns to the right, we calculate the total transmittance 
T = J2n=i l®i P- Physically, this corresponds to connect¬ 
ing the left and right edges of a finite network to uniform 
multi-mode waveguides, which act as scattering leads. 
The total transmittance is then measured as a function 
of the gain/loss parameter 7 (which is assumed to be ex¬ 
ternally tunable, e.g. by electrical or optical pumping). 

The variation of the transmittance with 7 is shown in 
Fig.[§a). Here, we take strip length M = 10 and width 
N = 10, and set the network links to (p = —0.257r, cor¬ 
responding to one of the Dirac points. The transmission 
is found to be peaked at certain values of 7; the peak 
positions depend on the choice of periodic or Dirichlet 
edge conditions, as well as the strip size. Fig. [^b) plots 
the values of 7 at the transmission peaks, versus the strip 
length M- For periodic edges, the peaks can be fitted to 


kch) 


(m -|- l/2)7r 
M 


TO G Zq , 


(15) 


where fee(7) is the tachyonic critical wavenumber given 
by Eq. (13). For Dirichlet edges, on the other hand, the 
transmission peaks can be fitted to 


kch) 



TO G Z+. (16) 




FIG. 5. (a) Plot of the transmittance T versus the gain/loss 
parameter g, for a strip of length M = 10 and width N = 
10, and quasi-energy p = —0.257r. Two cases are shown: 
periodic edges, and terminating (Dirichlet) edges, (b) Plot of 
the values of g at the transmittance peaks, versus the strip 
length M- The strip width is fixed at A = 10. Open circles 
are for a strip with periodic edges, and closed symbols are 
for a strip with Dirichlet edges. The curves show p red ictions 
made with the tachyonic dispersion relations, Eqs (15l-(16). 


The accuracy of these fits can be seen in Fig. [^b), by 
comparing the solid curves, which are produced from 
Eqs. (15) and (16), to the circles, which correspond to 
the numerically-obtained transmission peaks. 

The relationship between the transmission peaks and 
the critical wavenumber of the tachyon modes can be 
understood as follows. At the mid-gap quasi-energy 
(p = —0.257r, only the tachyon modes are propagat¬ 
ing; the other modes are evanescent, and thus incapable 
of forming standing-wave resonances. The transmission 
peaks occur when the strip length M equals (to -|- \)/2 
tachyon mode wavelengths, which allows for the largest 
intensity at the output column (intensity anti-node) rel¬ 
ative to the input column (intensity node). Although the 
tachyon modes have real wavenumbers, they do not over¬ 
lap exactly with the input amplitudes, so the presence 
of gain in the network results in overall amplification, 
i.e. transmission peaks higher than unity. For periodic 
edges, Eq. (15) follows from taking the tachyon modes 
to be plane waves propagating parallel to the strip. For 
Dirichlet boundary conditions, the tachyon modes must 
undergo reflections from the strip edges, and taking the 
lowest-order waveguide modes leads to Eq. (16). The 



























6 


dependence of this equation on both the width N and 
length M emphasizes the fact that the resonances arise 
from tachyon modes propagating in 2D, described by a 
PT-symmetric 2D Dirac equation. 

In conclusion, we have shown how an isolated tachy- 
onic dispersion can be realized in a photonic network, be 
it based on optical ring resonators or microwave trans¬ 


mission lines. A possible use of this tachyonic dispersion 
may be delay lines of wide tunability with low loss. 
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